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1.  General  Introduction 


Solid  vs.  Liquid  Triggering  in  Energetic/Explosive  Substances 

Triggering  (initiating)  is  one  of  the  key  processes  in  combustion,  explosion,  and  detonation  of 
energetic/explosive  materials.  By  controlling  the  triggering  mechanisms,  these  strong  processes 
can  be  made  more  probable  or  suppressed  if  necessary.  In  order  to  get  such  control,  one  has  to 
clearly  understand  the  physical  nature  of  the  triggering  mechanisms  in  different  substances. 
Depending  on  the  substance,  the  triggering  mechanisms  can  be  quite  different. 

The  processes  of  triggering  hot  spots  by  adiabatic  loading  of  liquid  energetic/explosive  materials 
have  been  explored  quite  well  (e.g.,  the  classical  monographs  in  references  1-4).  It  was 
established — both  theoretically  and  experimentally — that  the  key  role  in  these  processes  is 
played  by  various  gaseous  bubbles  (figure  1).  Because  the  bubbles  are  elastically  much  softer 
than  the  surrounding  energetic/explosive  liquids,  they  reach  much  higher  temperatures  than  the 
liquids  under  the  same  impact  pressure.  Thus,  they  are  the  hot  spots  for  initiating  corresponding 
exothermal  chemical  reactions. 

In  several  respects,  the  impact-induced  adiabatic  triggering  in  energetic/explosive  solids  differs 
from  the  corresponding  liquid  materials.  Three  essentially  new  circumstances  appear  when 
dealing  with  triggering  hot  spots  in  solid  energetic/explosive  substances: 

1 .  At  equilibrium,  liquid  substance  is  able  to  sustain  hydrostatic  loading  only,  whereas  all 
spectrum  of  nonhydrostatic  loading  becomes  possible  in  solids. 

2.  Also  at  equilibrium,  the  gaseous  bubbles  in  liquids  are  always  spherical  due  to  the  surface 
tension.  In  solids,  the  dissolved  gases  or  vapors  accumulate  within  the  pre-existed  defects 
of  various  shapes:  voids,  cracks,  intergrain  spaces,  etc. 
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3.  The  energetic/explosive  liquids  completely  transmit  the  externally  applied  pressure  to  the 
gaseous  bubbles  dissolved  within  them.  When  dealing  with  gaseous  bubbles  in 
energetic/explosive  solids,  the  solid  matrix  withstands  itself  to  the  externally  applied 
transmitted  loads. 

In  our  studies,  we  sequentially  analyze  each  of  these  factors  and  their  implications.  In  particular, 
we  dwell  on  a  detailed  exploration  of  the  effect  of  adiabatic  loading  of  the  solid  media  with  a 
spherical  hole  possessing  another  substance  within  itself.  This  study  gives  theoretical  and 
engineering  insight  into  the  mechanisms  of  generating  hot  spots  in  solid  energetic/explosive 
materials. 


2.  The  Exact  Nonlinear  Systems  of  Adiabatic  Loading  and  Its  Linearized 
Version 


2,1  Geometrical  Settings 

The  geometry  of  the  system  to  be  analyzed  is  shown  in  figure  2.  The  entire  space  is  referred  to 
as  the  Cartesian  coordinate  x' ,  with  the  basis  covariant  basis  x,  *  and  contravariant  basis  x  . 


Figure  2.  The  geometry  of  the  system. 


Although  the  systematic  use  of  “tensorial”  grammar  and  language  are  not  particularly  important  when  using  the  Cartesian 
coordinates,  the  formulas  look  much  more  elegant  and  well  structured  when  this  grammar  is  respected.  This  grammar  allows  the 
Cartesian  coordinate  system  to  be  replaced  with  others.  The  reader  can  ignore  all  the  tensorial  adjectives  without  losing  essential 
physics.  Note  that  in  the  Cartesian  coordinates  systems,  the  covariant  and  contravariant  components  of  all  tensors  are 
numerically  equal. 
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The  Lagrangian  (material)  coordinates  of  the  body’s  particles  in  the  unstressed  reference 
configuration  coincide  with  the  corresponding  coordinates  of  the  space.  We  use  the  notation  of 
r(x')  =  x'xi  for  the  radius  vector,  with  the  pole  at  the  origin  of  the  coordinate  system.  Let 
n  =  nx  =  n  x;.  be  the  components  of  the  unit  normal  to  the  hole’s  boundary  in  the  reference 
configuration.  For  the  spherical  hole  centered  at  the  origin,  the  components  of  the  unit  are  given 
by  the  formula  n‘  =  x‘  |x‘  (no  summation!).  Let  u(x)  =  u‘ (x)xi  =  w(.(x)x'  be  the  displacements  of 
the  material  point  with  the  Lagrangian  coordinates  x' . 

The  analysis  in  this  report  deals  with  thermoelastic  stresses  within  an  infinite  isotropic  matrix 
possessing  a  hole  and  exposed  to  the  specified  displacement  gradients  Ktj  at  infinity.  In  other 
words,  the  displacements  at  infinity  are  given  by  the  following  formula: 

it  — »  KjX1  at  |  xm  |  — »  oo  .  (1) 

2.2  Some  Necessary  Elements  of  Thermodynamics 

In  our  notation  of  thermomechanical  notions,  we  follow  the  monographs  in  references  5-8. 

Local  thermodynamic  state  of  any  elastic  solid  can  be  fixed  by  the  pairs  (ui  . ,  //)  or  (ui  jf0), 

where  ui  .  is  the  displacement  gradient,  77  is  the  specific*  entropy  density,  and  0  is  the  absolute 

temperature.  This  statement  means  that  all  other  thermodynamic  functions  are  uniquely 
determined  by  any  one  of  these  pairs.  It  is  particularly  true  regarding  the  internal  energy  density 
e  and  the  free  energy  density  1//  .  If  the  energy  density  e  is  given  as  a  function  of  (uj  . ,  77)  or  if 

the  energy  density  1//  is  known  as  function  of  (ui  j,9) ,  i.e., 


e  =  e(uu, 77),  y/  =  y/(utJ,0)’ 


(2) 


then  any  other  thennodynamic  function  can  be  calculated  only  by  differentiation  and  solving 
some  algebraic  equations. 

The  Piola-Kirchhoff  stress  tensor  pJ‘  can  be  treated  as  a  function  of  (ul  . ,  77)  -  pjl  (uk , ,  77)  or 
(uk , ,  0)  -  pjl  (uk  j ,  6)  .  Of  course,  it  is  silently  assumed  that  the  two  functions  are  completely 
different.  Sometimes,  in  order  to  avoid  possible  confusion,  we  use  the  notation  p  n  for  the  first 
function  and  gp]'  for  the  second  function. 


The  first  and  second  laws  of  thennodynamics  imply  the  following  identities: 


ij(  ,  Se(ukl,?j)  de(ukl,7]) 

nP  (ukl,p)  = - - - ,  0(ukJ,Tj)  = - - 


du. 


dij 


(3) 


and 


In  the  following,  the  adjective  “specific”  means  “per  unit  volume  in  the  reference  configuration.” 
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(4) 


y,  m  Sy/(ukl,0)  dy/(ukl,0) 

eP  (ukl,0)  = - - - ,  rj(ukJ,0)  = - 


du 


UJ 


de 


The  displacement  gradients  and  the  absolute  temperature  are  usually  the  most  convenient 
quantities  for  direct  experimental  measurements.  That  is  why  the  function  i//(ukl,0)  plays  such 
an  important  role  in  experimental  and  theoretical  thermodynamics.  Its  first  derivatives  can  then 
be  interpreted  with  the  help  of  the  thermodynamic  identities  in  equation  4.  It  is  much  harder  to 
measure  the  Piola-Kirchhoff  stress  tensor  and  the  entropy  density. 

The  second  derivatives  of  the  free  energy  define  the  so-called  instantaneous  isothermal  elastic 
modules  e  ciJkl ,  the  instantaneous  tensor  of  thermal  expansions  j3iJ ,  and  the  instantaneous  heat 
capacity  cu  at  fixed  displacements: 

Jjki  _  d' ,  0)  oj  _  d2y/{umn ,  0)  dp"  (um  n ,  0) 

]C  —  _  _  ,  p  — 


duijdk,i 


du,  £ 


d6 


Qd2V(um,n£)  drj(u  ,0) 
c„  =  -0 - - - =  0- 


d0l 


d0 


(5) 


“Instantaneous”  means  that  these  quantities  change  when  the  thermodynamic  state  changes  itself. 
The  traditional  constant  tensor  of  the  linear  elasticity  is  just  the  values  of  these  tensors  in 
reference  (usually  stress-free  or  hydrostatic)  state. 

The  internal  energy  density  function  e(uk  l ,  rj )  plays  a  particularly  important  role  when  dealing 

with  adiabatic  processes  in  which  the  entropy  density  remains  unchanged  in  each  material 
particle.  Its  first  derivatives  are  given  by  the  thennodynamic  identities  in  equation  3. 

The  instantaneous  adiabatic  elastic  modules  „  clikl  are  defined  by  the  second  derivatives  of  the 
specific  internal  energy  density:  nc  ,ki  =  d2e(um  n,r/)/ dui  duk  t . 


The  first  and  second  laws  of  thennodynamics  imply  certain  thermodynamic  identities  binding  the 
second  derivatives  of  the  two  main  thermodynamic  potentials.  In  particular,  these  identities 
imply  the  following  useful  thermodynamic  relationships: 

ciJkl  =ecijk! +  —/3ijj3kl  ,  (6) 

c„ 

and 

d2e(uu,rj)  _  0  d2e(ukJ,ij)  _  0 

dr/2  cu  ’  dr/dui  j  cu 
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2.3  The  Exact  Nonlinear  System  of  Equilibrium  Equations  for  Isentropic  Loading 

We  assume  that  the  original  stress-free  state  of  the  body  has  a  uniformly  distributed  entropy 
rj  =  r/0 .  We  will  also  assume  that  in  the  loading  processes,  the  entropy  distribution  remains 
unchanged.  The  following  equilibrium  conditions  should  then  be  satisfied  within  the  bulk  of  the 
body  and  at  its  stress-free  boundary,  respectively: 


d  Se(umn,T]0) 
dx 1  dui  j 


(8) 


and 


Se(um„,rio) 


du 


i,j 


rij  =  0 . 


(9) 


After  solving  equations  8  and  9,  i.e.,  after  determining  the  vector  field  of  displacements  ut(x)  ,the 
equilibrium  distributions  of  the  stress  tensor  p  j‘  and  the  absolute  temperature  can  be  calculated 
with  the  help  of  the  thermodynamic  identities  in  equation  3: 


Se(ukJ,rj0) 

duu 


(10) 


and 


ad@(llk,l)  ~ 


ge(“t,/.7o) 

dr/ 


(ID 


3.  The  Equilibrium  System  in  Linear  Approximation 


Assume  that  the  loading  experiment  can  be  adequately  described  in  the  framework  of  quadratic 
approximation  of  thermodynamic  potentials  near  the  stress-free  configuration.  In  particular,  this 
means  that  the  free  energy  density  can  be  approximated  by  the  following  relationship: 


1  dV(O,0o)  ,  SV(O,0o)  1  52y/(umn,e0) 


duijd. 

iuk  1  ~  ui 

j  'J  ’  dutJd0 

rijklu 

0^  Ui, 

iuk,  +  BiJuilT--^T2 

J  lJ  2  e0 

dO2 


(12) 
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where  tensors  0C,jkl  =  0cljkl (0,0Q)  ,  Bkl  =  /?w(0,#0) ,  and  Cu  =  cu  ( 0, 0O )  play  the  roles  of  the 
second-order  instantaneous  tensors  in  the  ground  configuration;  T  =  0-6o  is  the  notation  of 
deviation  of  the  absolute  temperature. 

Equation  12  leads  to  the  following  linear  expansion  of  the  stress  tensor  p" : 

P,\K,„T)  =  ,C’“utJ+B’T,  (13) 

the  bulk  equilibrium  equations, 

«£*%,+«%-<>,  (14) 

and  boundary  conditions  at  any  traction-free  interface  S  , 

(.CVu^+tfiT^nj-O.  (15) 


When  dealing  with  adiabatic  loading,  it  is  more  convenient  to  work  with  the  following  bulk  and 
boundary  equations  of  equilibrium: 


v 


=  0, 


(16) 


and 


n ,  =  0 . 


(17) 


The  corresponding  stress  tensor  and  temperature  change  are  given  by  the  following  formulas: 


ad  l 


Cijklu 

Uk,l 


(18) 


and 


de(0,?i0)  ^ 

ad 1  ~  Uk,l 

dmk. 


(19) 


With  the  help  of  the  identities  in  equation  7,  equation  19  can  be  rewritten  as 


ajT  ~—Bklukl . 
c. 


(20) 


Using  relationships  equations  6  and  20,  equation  18  can  be  rewritten  as  follows: 


C»“+k!LB"B 

c 


uk,i  =  oCiJk\,  +  jr# , 


(21) 


as  it  follows  from  equation  13. 
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Isotropic  Elastic  Substance 

In  the  case  of  isotropic  elastic  substance,  the  tensors  eC'jkl  and  Bij ,  by  definition,  have  the 
following  particularly  simple  structures: 

e  CiJkl  =  A .g5ij5kl  +  n0  (S‘kSn  +  S“SJk ),  BlJ  =  -K0a5" ,  (22) 

where  Ae  and  fie  are  the  Lame  isothermal  elastic  modules,  Ke  =  Ae  +  2 //^  /  3  is  the  isothennal 
module  of  volumetric  compression,  and  a  is  the  coefficient  of  thermal  expansion.  The 
definitions  in  equation  22  imply 

p,J  =  A9Sijukk5kl  +  2 jugu^  -  KgadijT .  (23) 


Some  formulas  of  thermodynamics  of  isotropic  solids  become  more  convenient  when  using  the 
isothermal  and  isentropic  Poisson  ratios  ve  and  v : 


/L 


JL 


ve  =' 


2 (K+Me)’  ’  2(A,+ft,)' 

Formulas  in  equation  24  imply  the  following  useful  identities: 

2//  1  +  Vg 


,  _  ..  2ve  ,  ..  2vn 

^8  Be  ,  ^  Ml 


1-2  v,  "  -"l~2vn 


,K 8  = 


- 

3  l-2v.  3  l-2v. 


K0=  xs  +  U,  =  i±r«L,  *  =  4  + 

3  3  l-2v„  ’  7  3  n  3  1  - 2v 


(24) 


(25) 


Using  equation  22,  we  can  rewrite  equations  20  and  21  as  follows: 

.J  =  ~Ksau‘t, 


c 


and 


"  W  +  2/W  j  -  ^  ^  • 

With  the  help  of  equation  22,  equations  6  and  7  can  be  rewritten  as  follows: 


(26) 


(27) 


=  0Cijkl  +  ^-K2a2SijSkl  =  Ad5i]5kl  +  p0  ^ 


0, 


A9  +—K29a2 

\  Cu  j 


SijSkI  +  jUg(SikSjl  +SilSjk), 


(28) 
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and 


d2e(0,7o)  _  52g(0,7o)  _  #o  ^  aSv 

dr  Cu  ’  5//5u;  /  C„ 

Introducing  the  isentropic  Lame  modules  /l;?  and  according  to  the  formulas 

Q 

—  k 0  "t-  —  K-q  cc  and  1 1 ^ i 9 , 


(29) 


(30) 


we  can  rewrite  equation  28  as  follows: 


n  Cijkl  =\5ij5kl  +  r  ( 8ik8il  +  S"  Sjk ) .  (31) 

According  to  equation  27,  the  stresses  in  the  adiabatic  loading  take  the  familiar  form  as  follows: 

ay  «\SiJukkSk'  +  2^''.  (32) 


Equation  30  implies 


K„  =K9  +  ^Kla\  —  =  C-‘  +  K°6,/r 
C ..  9  Ka  C 


(33) 


Also,  we  would  like  to  indicate  the  thennodynamic  identities  which  include  the  specific  heat 
capacity  Cp  at  fixed  pressure: 

c,  -  c,  =  KAa\K,  =Ke+  A Ky ,  A  s  A .  (34) 

Kg  L1( 


4.  Temperature  Distribution  in  Isotropic  Elastic  Space 


4.1  Adiabatic  Nonhydrostatic  Loading  of  Unbounded  Isotropic  Space 

This  case  is  based  on  the  analysis  of  equation  27,  which  can  be  rewritten  in  the  following  form: 

+2/y<<'i>  =,y+KA,TS'.  (35) 


Equation  35  implies  the  following: 


3K, 


k. 

'  adP  .k 


+  aadT  = 


3  K, 


k. 

'  adP  .k 


(36) 
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Combining  equations  26  and  36,  we  get  the  following: 

_  OqOC  ad  P  .k  _  QqOC  ad  P  .1 

ad  " ~Cu+60a2Ke  3  "c7~' 

Using  equations  34  and  37,  we  have  more  convenient  fonnulas: 

T  _  Oj/y-  Kg  adp\k  _  60a  1  +  vg  1-2^  adpkk 
ad  "  C„  Kn  3  '  C„  l-2v,  1  +  v,  3 


(37) 


(38) 


4.2  Formal  Statement  of  the  Boundary  Value  Problem 

Summarizing  the  previously  mentioned  relationships,  we  arrive  at  the  boundary  value  problem 
(BVP)  in  linear  approximation. 

The  bulk  equations  of  equilibrium: 

~\8iJ8u  +  nn  ( SikSJl  +  8“8jk )]  ukjl  =  0 .  (39) 

The  condition  at  the  hole’s  boundary: 


' \8j8kl  +  nn  [8,k8il  +  8“8Jk )' 


*k,i 


rij  =  0 . 


(40) 


The  conditions  at  infinity: 


ut  ->  KjXJ  at  |  xm  |  — >  oo . 


(41) 


After  solving  the  BVP  in  equations  39-41,  the  impact-induced  temperature  change  T  can  be 
found  with  the  help  of  equation  26. 

4.3  Formal  Solution  of  the  BVP  in  Equations  39-41 

We  will  be  looking  for  a  solution  of  the  system  in  equations  39-41  in  the  following  solution: 


ui  =  KaxJ  ~  4(1  -  vrj  )Dip 


rn 

vu 


r 


+  D 


pq 


l 


r  +  C- 
rj 


1 


,iPq  =  ^+A(\-v^Dip  -x* 


+  . 


\(8pqxi  +8lpxq  +8‘‘xp) 

(3 

C— _ — 

^5  3 

—  xp  xq  xi 

fc— )1 

^7  5  f 

{  J 

L  r  r  ) 

l  r  r  )\ 

(42) 


where  C  and  D j)q  are  the  coordinate-independent  scalar  and  symmetric  tensor,  respectively,  and 
r  =  ^8ijxixJ  and  v  =  XJ 2(\  +  //;?)  are  the  adiabatic  Poisson  ratios. 


A  somewhat  routine  but  cumbersome  calculation  leads  to  the  following  formulas: 
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u,=KttxJ  +4(1  -vn)D- 


+  D. 


pq 


(3  1  ) 

(  15  3  ')] 

\[Spqxi+Spxq  +8qxp) 

C—~  — 

|  -8lrxpxqxr\ 

r  —  i 

^  7  c  f 

V  1  1  1  j 

v  r  >'3  J 

l  r  r  J  J 

(43) 


and 


,,,=^-4(1 

(*%  +8*8)  +Sfddc^-^\  +  3d.xr(Sp‘‘dimx-  +S?x«  +8?x’)\±f-Cy\ 
-3  (S’x'SS  +x*5]5,x  +xPtfSv){c^-^-l5x>x°8^8j9tf[±-C^\ 


+D  { 

pq  ' 


(44) 


The  formula  of  the  stresses  is  even  more  cumbersome: 


2v„ 


1 

—Pj  =  ir~ 
M  1-2  v„ 


-8„ 


4 +2(1 -2  v,)D*^ 

r 


■12  vS^tfx'x"- 
r 


+  2Km-4(l-vn) 


(  3  1  ^ 

Dik  -r+  j 

1 

. 

(  |  CO 

cf 

1 

L  v  r  J 

7  J 

-4(1  -vn)Djk-Simxkxm 
r 


+  2D 


P‘1 


(<>  "4 + 24  '■<>;' )  |c + - +J  -  34,,-t"  ( <>  "4-'J  +24V) 


,7  5  , 

V  r  r  j 


+  2D., 


1 544.x  Wx*|c2r-- L 


9  7  , 

r  r  j 


-3(2  S(pxq)Simxm  +  xpxq8i y) 


'cA_lA 

,7  5  . 

V  r  r  j 


(45) 


Straightforward  calculations,  based  on  equation  43,  show  that  the  asymptotic  condition  in 
equation  41  at  infinity  is  satisfied  automatically.  More  cumbersome  calculations  show  that  the 
bulk  equilibrium  in  equation  39  will  be  satisfied  automatically  as  well.  Lastly,  straightforward 
calculations  with  equation  45  show  that  the  boundary  condition  in  equation  40  will  be  satisfied  if 
the  constant  C  assumes  the  value  of 


C  =  -R2 
5 


(46) 


and  the  symmetric  tensor  satisfies  the  following  linear  algebraic  system  of  equations: 


^2-5v  l-5v  ,  ^ 

- n-8pq8..  + - q-28(p8q) 

5  v  5 


J 


n  Ri  °° 

D  =  —  p. 

2/u,J 


pq 


(47) 
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where  pf  is  the  value  of  the  stress  tensor  at  infinity 


00 

Pa  =  2P 


l-2v 


~KA  +  Km 


(48) 


Equation  48  implies  the  following  useful  relationship: 


1  -  2v„ 


2//  l  +  v 


(49) 


Equation  47  has  the  following  solution: 


Dm=~ 


R 


2-5v„ 


»  4(4-5vJ(7-5^) 


00/C  C* 

PkSpc  + 


fc 


M  4(7-5 v,) 


/Ur 


(50) 


4.4  Temperature  Distribution  Around  a  Spherical  Hole 

Combining  the  formula  of  the  displacement  gradient  in  equation  44  with  equations  46  and  49, 
we  get 


*0  =<  +2(1 -2^)^ -3 


if  3  A 

Sqp-^xpxq 
\  r~  j 


(51) 


Combining  equations  49  and  51,  we  arrive  at  the  following  formula  of  the  local  volume  change: 


2//  \  +  v  r 


\ 


Spq-^xpxq 
v  r2  j 


(52) 


Using  equation  50,  we  can  rewrite  the  formula  of  the  divergence  as  follows: 


k 

U..u  = 


1  l_2v  x  R3  1  5(1  -2v) 

P.k  +■ 


3 

'/'I  xU  _  J  Jvi 


2//  l  +  v 


- xl  X 

M  2{7-5vn)[d  rXX  , 


(53) 


Finally,  substituting  equation  53  in  equation  26  and  using  equation  25,  we  arrive  at  the  following 
formula  of  the  temperature  change: 


ad 


T  =  _a0o  l  +  v0  l-2v„ 


3CU1-2V0  \  +  vn 


SiJ  +  5R  l  +  Vq  f 


r3  7  -  5v„ 


3  .  A 

S,J-^x'xJ 
r  j 


Pa- 


(54) 


With  the  help  of  identities  in  equations  25  and  34,  equation  54  can  also  be  rewritten  as  follows: 
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(55) 


ad 


T  =  - 


a^0  Kq 


3 C,  K 


u  rj 


r3  1  +  vn  f 

Su+5 


r  7  -  5v„ 


3  .7 

S]  -4rXlXJ 
v  r  J 


a  6,, 


3  C 


£"+5 


/?3  l  +  v„ 


A 


r  7  -  5v„ 


y 


v  r  y 


Pi!  • 


This  is  the  main  result  of  this  report,  which  is  discussed  in  the  next  section. 


5.  Discussion  and  Conclusion 


Equations  54  and  55  are  quite  instructive  and  transparent.  The  first  terms  give  the  spatially 
uniform  adiabatic  temperature  increase  Tuniform  that  would  exist  in  the  absence  of  the  hole  in  the 

infinite  elastic  media: 


T 

uniform 


a6. 


3  C 


°  Pa  = 


a9o  Kg 
3  C„  Kn 


Pa 


(56) 


The  second  term  gives  a  spatially  nonuniform  additional  temperature  increase  Tlocal  due  to  the 
presence  of  the  hole: 


T 

ad  local 


=  o  1  +  4  R 


3  l-5vv  r 


v  q  y 

Sy-ArXiXj 


V 


) 


(57) 


The  temperature  field  Tlocal  is  localized  in  a  small  vicinity  of  the  hole  and  decays  as  the  inverse 

cube  of  r  distance  from  the  center  of  the  hole.  However,  at  the  hole’s  boundary  at  r  =  R,  both 
temperature  fields  TV,  and  Tlocal  are  the  same  order  of  magnitude. 


Remarkably,  the  local  temperature  field  Tlocal  identically  vanishes  when  the  solid  undergoes  the 
pressure-like  external  loading  at  infinity,  i.e.,  when  p'j  =  -pextSj  .  On  the  other  hand,  the 
uniform  temperature  field  Tunjform  vanishes  at  shear-like  loadings,  i.e.,  when  pf  =  0 .  In  this 
case,  the  local  field  Tlocal  not  only  dominates  but  becomes  the  only  existing  temperature  field 
induced  by  the  adiabatic  impact.  In  this  case,  it  can  be  rewritten  in  the  simpler  form: 


T 

ad  local 


ag,5(n-v,)^ 

Cr  7~5vn  r>  P‘ 


(58) 
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where  r‘(xk)  =  x  /  r  are  the  Cartesian  components  of  a  unit  vector  having  the  same  direction  as 
the  radius-vector  f{xk )  drawn  from  the  center  of  the  hole  to  the  point  x  . 


Equations  57  and  58  become  even  more  transparent  in  terms  of  the  three  principal  stresses,  pL 
of  the  tensor  p7 ,  corresponding  to  the  three  principal  directions  with  the  unit  vectors  l]  .  In 
terms  of  the  principal  quantities,  the  tensor  p..  can  be  presented  in  the  following  form: 

pt  =  ipV5-  (») 

L= 1 


Introducing  the  directional  cosines  cos  0L  =  rl' ,  we  can  rewrite  equation  57  as  follows: 


T,ocal(X)  = 


a0o5(i  +  vn)R 
C„  7  -  5v„ 


3  3 


-'£pLsm20L  . 


L= 1 


(60) 


In  particular,  for  the  uniaxial  stressing  along  the  axis  x 1 ,  equations  56-60  give  the  following 
formula  of  the  spatial  distribution  of  the  impact-induced  temperature  field: 

(  n3  A 


ad 


T=-  5a0° 


1  +  V 


7-5^ 


1  -  ^  sin 2  6X 


(61) 


On  the  boundary  of  the  hole,  equation  61  gives  the  following  formula: 


5(l  +  vn) 

CP 


p'  cos2  . 


(62) 


All  these  remarkable  qualitative  facts  should  be  taken  into  account  when  discussing  generation 
of  hot  spots  in  solids  under  the  action  of  adiabatic  impacts.  The  explicit  solution  can  also  be 
used  for  verification  of  numerical  codes  and  planning  experiments.  At  this  stage,  it  is  important 
to  merge  theoretical  modeling  with  computer-based  modeling  and  physical  experiment. 
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